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<N . Abstract 

& - 

Following the general results on the relationships about Backlund transforma- 
tions (BTs) and exact discretisation given in a previous work [12], we consider 
[ — . ' the Ablowitz-Ladik hierarchy and a corresponding family of BTs. After dis- 

cussing the boundary conditions, we show how to get explicit transformations. 
The Hamiltonian properties of the maps and of the discrete flows are examined. 
The conditions on the parameters of the map giving exact discretisations are 
discussed. Finally, analytical and numerical examples are given. 
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on . 1 Introduction 
o 

m 

There exist a fairly large number of papers explicitly or implicitly related to the subject 
of Backlund transformations (BTs) applied to the Ablowitz-Ladik (AL) model [TJ[2]: 
here we just recall the paper by Suris [9J and the more recent one by Tsuchida [ID] . 

These two authors adopted the original approach of Ablowitz and Ladik [3], based 
on the zero-curvature representation of the model. More explicitly, suppose that our 
lattice model possess a Lax pair given by (Lfc(A), Mfe(A)), where k is the label for the 
lattice site and A is the spectral parameter, i.e. an arbitrary constant independent of 
k and t. Then, the linear problem associated to the model is given by [2]: 

£*(A)Vfc(A) = ^+i(A), = M fe (A)^(A), (1) 

To the solution ^(A) of ([I]) we can associate another solution ip^ of the same 
problem multiplying ipk by a suitable matrix Dk, usually called Darboux or dressing 
matrix. The set of equations obtained [3J, that is 

£*(A)V*(A) = ^+i(A), MX) = D fe (A)V fe (A), (2) 
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are the natural time discretisation of the system ([I]) if we consider the new solution 
iftkW as the old solution ^(A) but computed at the next time step. The compatibility 
condition implied by <^ is: 

L k (X)D k (X) = D k+1 (X)L k (X). (3) 

The above procedure is broadly accepted to be a good starting point to obtain 
discretisations, but from this point on, there is a bifurcation in the literature results: 
there are various expressions for the dressing matrix D k , in particular for its depen- 
dence on the spectral parameter, and the solutions of the over determined system 
are presented in different forms. In [9] Suris showed how the non local maps found by 
Ablowitz and Ladik [3] can be factored in the composition of simpler and local trans- 
formations, although defined only implicitly. Then Tsuchida, in his interesting paper 
[TO] , considers the relation [3] for a given lattice system, not necessarily given by the 
AL model. Also he considers the expression obtained taking the determinant on both 
sides of[3j 

det(L fc (A)) det(£> fc (A)) = det (£>*+! (A)) det(L fc (A)), 

showing how this is sufficient, under certain conditions, to obtain local expressions for 
the maps representing the discretisation of the underlying system. The author is aware 
of the fact that the discretisation obtained with this method are indeed BTs. However, 
he seems to not recognize that the algebraic equations solved by certain particular 
combinations of the dynamical variables appearing in the matrix D k (called "auxiliary 
variables" in [ID]), can be found by evaluating the determinant of the dressing matrix 
at a certain constant value of A and then setting the result equal to zero. 

The relevance of this fact comes from the observation that it allows to obtain 
explicit transformations, as shown for example in [7] , [8] , [UJ [12] . This may imply an 
improvement in the computations of the discretisation scheme. 

Basically the maps of Ablowitz and Ladik are non-local and implicit, the maps of 
Suris and Tsuchida are local and implicit, here we show how to get explicit maps, in 
general non local. Further, we show how our formulation gives a remarkable Hamil- 
tonian characterization of the discrete flow. Indeed, since the BTs usually involve the 
whole hierarchy of commuting flows, it is reasonable to ask which of the flows of the 
hierarchy we are really discretising. More importantly, the question could be about 
the possibility to arrange the scheme so to discretise only a particular flow, e.g. the 
one corresponding to the physical Hamiltonian of the model. 

We will try to answer these questions following the results of [12], where has been 
shown that the BTs for an integrable systems are deeply connected to the Hamiltonian 
structure of the system itself. Indeed they represent the integral curves of another in- 
tegrable non-autonomous system that shares the conserved quantities with the original 
model: this observation can be used to obtain an explicit formula for the Hamiltonian 
interpolating the discrete flow defined by the iterations of the BTs. Furthermore, for 
any fixed orbit, the discretisation lies on a suitable (linear) combination of the integrals 
of motion: by choosing in a proper way the parameters appearing in the transforma- 
tions it is then possible to make the discrete flow and the continuous one parallel for 
all time. 
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The paper is organized as follows: in section [2] we briefly recall the results from |12j ; 
in section [3] we describe how to properly define the overall Lax matrix under periodic 
boundary conditions or in the case of an infinite support; in section H] we present the 
core findings of the work, that is the formulae for the explicit discretisation of the 
model and finally, in section we discuss some analytical and numerical experiments. 



2 Hamiltonian flow associated to BTs 

First of all, let us collect some known facts about BTs. Suppose to have an integrable 
system with N degrees of freedom and N quantities {Hi\f =1 conserved and in involu- 
tion. Denoting by (pi,ft)£Li the canonical coordinates, the Poisson structure is given 
by: 

{Pi, Qj} = ftj, {Pi,Pj} = {ft, Qj} = 0. 



A diffeomorphism from the phase space to itself, say 



(Pi,qi)=> (Pi,Qi), where { r * (4) 

v ft = 9i[p,q), 

is a BTs if the following two conditions are satisfied: 

• The transformations are canonical: 

{pi, Qj} = {PiiPj} = {ft, Qj} = 0. 

• The functions of the phase space given by the conserved quantities are invariant 
under the action of the map: 

Hi(p,q) = Hi(p,q). 

To simplify the notation we will drop hereafter the subscripts in the independent 
variables as above. If the transformations are parametric the previous definition re- 
mains the same. In the rest of this section we assume to have parametric transforma- 
tions with a parameter /x; further, we assume that the BTs are connected to the identity, 
so that (ft, q)\ = (p, q). As shown in [12], these assumptions are not restrictive. 

One of the most interesting features of parametric BTs is the spectrality property 
[7] : it is deeply related to the Hamiltonian structure of the transformations themselves. 
Let us call the parameter of the transformations fi. Since the maps are canonical for 
every value of the parameter /i, there exists a generating function F such that [S]: 

dF(q,q, /i) 
~ _ dF(q,q,fi) 

Pi o~ 



Pi = fi(p,q), 
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The function F is a function of the 2iV variables q^fLi and of the parameter ji. Let 
us call G(q, q, fi) the partial derivative of F with respect to fi. We say that the BTs 



through the relations fll]), that is 



possess the spectrality property if, expressing G as a function of the variables (pi, qij{ : 



$(P,9,A*) = G(q,g(p,q),fi), 
the function $ commutes with all the conserved quantities of the system: 

{<S>(p,q, f i),H i {p,q)} = 0, i = l,...,N. 

An immediate consequence of this property is that the BTs are the integral curves, 
whit parameter /i, of the non autonomous Hamiltonian system governed by the function 
q,fi) [12], that is 

dq~i = d$(p,g,/x) 

9fi dpi ' _ 

dpi _ d$(p,q,fi) 
d/i dqt 

As a remark, we stress that q, fx), as a function on the phase space, is a function 
of the conserved quantities only [T2] : 

$(p,q,fj) = 

This result has important consequences as regards the computational application of 
the BTs. Indeed, denoting by and g| n ^ the n-th iteration of the maps, the discrete 
trajectories defined by the set {p^, q^}^ =l lie, for every value of T, exactly on the 
integral curves of the autonomous system 

dqt = &H{p,q,n) 

dt dpi ' , g s 

dpi _ dH(p, g, /x) 
dt dqi 

where the interpolating Hamiltonian H [12] is given by 

n(p,q,»)= r$(p,q,\)d\. (7) 
Jo 

The concrete possibility to use the BTs in numerical applications comes from the 
following observation: for any function $(p, q) on the phase space we can write 

Z={U,Z} = y j c k {H k ,Z}, where c k = -|- T Q(H, X)dX. (8) 

The functions c k are constants of motion, so they are constants along the trajectories of 
the model. Having in mind a numerical experiment, we must fix the numerical values of 
the initial conditions and those of the parameters. Then, the above formula tell us that 
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we can treat the parameters Ck just as constants: for any fixed set of initial conditions, 
the discrete trajectories coincide with the continuous ones given by the Hamiltonian 

N 

y / J CkH k , 

k=l 

where the numbers c k must be calculated according to the formulae (El)- 

Suppose now to have a multi-parametric set of BTs, depending both on the param- 
eter \x and on the set {Afc}| =1 , for some integer s > 1. The interpolating Hamiltonian 
([7]), as well as the constants Ck, will depend on this set of parameters. These additional 
degrees of freedom can be used to force the values of the Ck to be equal to some pre- 
assigned set. We will see how this method works in the next sections, where we will 
construct a two-parameter family of BTs for the AL hierarchy. 



3 The Lax matrix and a comment on the boundary 
conditions 

The AL model is described by the following equations of motion [2]: 

Qk = Qk+i + Qk-i ~ %qk ~ qkr k (qk+i + <?fc-i), , Q s 
f k = -r k+ i - r k _i + 2r k + q k r k {rk+i + 

where qk(t) and rk(t) are dynamical variables on a lattice. It is possible to consider 
both periodic boundary conditions, that is q k = qk+N and = r k +N for some iV G Z, 
and an infinite chain. For the infinite lattice we suppose that the variables possess a 
finite L 1 norm: 

^2 \Qk\ < oo ; ^2 l rfc ' < 00 

k k 

Let us discuss the definition of the Lax matrix of the model separately. 

Periodic case. In the periodic case, the whole hierarchy of the commuting con- 
served quantities can be found from the invariants of the Lax matrix L(\). This matrix 
is defined by the product 



L(X)=l[Lk(X), with L k (X)=^ k A ff iJ. (10) 



The determinant of L(X), given by HfeLi(l — qk r k), is a conserved quantity. The 
other N — 1 conserved quantities appear in the Laurent expansion of the trace of L(X): 

N 

Tr(L(A)) = H ^ N ' 2t > H = H N = 1. (11) 

i=0 
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The involutivity of these conserved quantities follows From the r(A, ^-matrix struc- 
ture satisfied by the Lax matrix: 



{L{X) f L(rj)} = [r, L(X) <g> L(rj)], 
where r(X,r]) is defined by [9]: 
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The above relations are equivalent to the following Poisson brackets for the dynam- 
ical variables of the model: 



{ik, Tj} = (1 - q k r k )5 kj , {q k , qj) = {r k , r,-} = 0. 
The flow corresponds to the Hamiltonian h given by: 



(13) 



k=l 



h = -H! - H N -i-2\og(det(L)), where Hi = y^g fc r fc _i, H N „ X = ^r fc g fc _i (14) 

k=l k 

Infinite chain. In this case the Lax matrix is given by the formula 

n 

l ^ = n 



(15) 



m-t-oo k=m 



whenever the limits exist. 

It is important to notice that the equation ([3]), defining the BTs, is homogeneous 
in L k : if we multiply the matrices L k by any scalar number a, the BTs are not affected 
by this operation. Let us denote by Ljj(A), i, j — 1,2, the elements of the Lax matrix 
(I15p . From the homogeneity of ([3]) it follows that the BTs are homogeneous function 
of degree zero of the elements Lij(X); for example for q k we will see that 



q n = q n (q,r, (pk{L it j(X k )),X k ),f n = r n (q,r, <p k (L itj (\ k )),\ k ) 



(16) 



where are a set of free parameters and (p k (Lij(X k )) stands for a set of functions of 
the matrix elements L it j evaluated in A = X k . The functions if k (L i j(X k )) must be all 
homogeneous of degree zero in the variables L it j(X k ), that is 



ip k (aL i} j(X k )) = ip k (Lij(X k )), 



Va e C. 



The Lax matrix (TTo*]) is well defined for |A| = 1 if the sequences {r k } and {q k } have 
finite £ l norm (see|A]). However, the functions ip k (Lij(X k )) are well defined even for 
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I Afc| 7^ 1 in all the points where they are continuous. To clarify this point let us consider 
the sequence of matrices T^: 



T N =f[ L k (X). 

k=-N 

The corresponding elements can be written as 

TV TV 
iN _ ST^ j(k) x 2k+l rpN _ ST^ ;(fc)\2fc 



rpN /( fc )\2fc+l rpN _ V /( fc )\ 2 

fc=-AT fc=-TV 
AT AT 

rpN _ /( fc ) \ 2fc _ /(fc) \2fc-l 

1 12 — / j t 12 A J J 22 — / l 22 A i 



(17) 



-TV fc=-TV 



where the coefficients ^ are functions of the dynamical variables. Consider, for exam- 
ple for |A| < 1, a homogeneous function y?(Z^ (A)) of degree zero. Due to the absolute 
convergence of for |A| = 1 and using the continuity of the function ip, we can write 

lim p(2#(A)) = lim ^(A^ +1 lf (A)) = <p( lim \ 2N+1 T^{\)) = ^(A)). 

TV— »oo J TV— >oo J TV— s>oo J 

The convergence of (p(T^(\)) for |A| > 1 can be proved similarly. In the case of an 
infinite chain, the BTs (|T6l) are then well defined also for values of the parameters A^ 
inside or outside the unit circle. 

In the following we will freely talk about the matrix L(X): the results will apply 
independently of the boundary conditions assumed. 

4 Discretisation of the Ablowitz-Ladik hierarchy 

Looking at the A dependence of L k (X) ffTU]) . we seek a dressing matrix D k of the form: 

D k (\)=( w ~ l{x+x ~ W) w a as) 

v ; V c k w{X^ + Xd k ) J v ; 

Here w is a constant, independent on the lattice label k, introduced for later conve- 
nience. 

As noticed in the introduction, to obtain explicit transformations it is important 
to have a dressing matrix that is singular when A assumes certain constant values. 
We will take a matrix D k (\) possessing two degenerate values of A, say Ai and A 2 : 
these values will appear also as free parameters of the transformations. Note that the 
constraints det(-Dfe(Ai)) = and det(D fc (A 2 )) = leave only two free variables among 
the four appearing in [18] We can parametrize D k as follows: 

A , AiA 2 (afcA2-/3fcAi) \\— \\ 



n f\\ — I w M>A(A 2 /3 fc -Aia fc ) A 2 /3 fc -Aia fc / 1q n 



AiA 2 (A 2 /3 fe -Aia fc ) V A Ai A 2 (A 2 /3 fc -Ai» fe ) 
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where a k and f3 k are two new variables. The matrices Dk(Xi) and D k (X 2 ) possess a 
kernel by construction, respectively given by 



By applying \Vt\ > to the relation (j3]) evaluated at A = Ai we obtain 

Lk{Xi)D k (Xi)\Ql >= = D k+1 (X x ) (^fe(Ai)l^ >) . 

Up to an overall factor, \Q\ > is unique. Then, from the previous equation, we can 
write 

L fc (Ai)|nj[ >=g k \nl +1 >, 

for some function g k . More explicitly, after eliminating gk, we obtain the recursion 

w(X 1 wr k + a k ) w w 2 (l - q k r k ) 
Ai(Aiw + q k a k ) Xiq k q k {Xiw + q k a k ) 

Exactly for the same reason, the variables f3 k must obey the recursion: 

^ ^ _ w(X 2 wr k + (3 k ) _ w _ w 2 (l - q k r k ) 
X 2 (X 2 w + q k /3 k ) X 2 qk qk(X 2 w + q k (3 k ) 

Remark. If we wish to construct a map preserving the conserved quantities, the 
boundary values for (3 k and a k are not arbitrary. 

Take for example the periodic case. The Lax matrix of the model is given by 
L = Ln...L\. Imposing a periodicity also on the dressing matrix, so that D^ +k = D k 
for every k, form the relations ([3]) it follows that 

L(X)D 1 (X) = D 1 (X)L(X). (22) 

Since the eigenvalues of £(A) are time- independent, the isospectral equation ( 122]) im- 
plies the preservation of the conserved quantities under the action of the map. The 
periodicity of D k means that the variables a k and f3 k are periodic as well: the relations 
(1201) and (I2TI) are then periodic continued fraction. This means that the variables a k 
and /3 k are solutions of quadratic equations. To find out these quadratic equations, 
we must recall that Di(Ai) and Di(A2) are degenerate matrices. Again, applying the 
kernels \Q\ > and > to the equation ( )22|) evaluated respectively at A = Ai and 
A = A2 we obtain 

LMACAiM >= = Di(Ai) (L(Xi)\Q\ >) , z = 1,2. 

This time we conclude that \Q\ > is an eigenvalue of L(Xi), so we can write 

L(A i )|fii>=7 i |nj> ) i = 1,2, 

from which, eliminating j i} we obtain the quadratic equations solved by a.\ and /3i 

Li, 2 (Xi)al + (Li,i(Ai) - L 2 , 2 (Xi))wa 1 - L 2)l (Xi)w 2 = 0, 
^i, 2 (A 2 )/3 1 2 + (A 2 ) - L 2i2 {X 2 ))w^ - L 2A (X 2 )w 2 = 0. 
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Vice- versa, if we start from these boundary values we will obtain a periodic dressing 
matrix, since the values a.N+k and Pn+u obtained from ( 120]) and (|2T|) are equal to those 
of a k and j3 k . Similarly, in the case of an infinite period, to take the limiting values of 
a k and f3 kl for k — > — oo, as the solutions of the equations (12"5|) . implies the equivalence 

lim£> fc (A)= lim D k {\), 

k— >oo k— >— oo 

which, in turns, implies the preservation of the conserved quantities under the action 
of the BTs. These results agree with the observations made by Tsuchida about the 
boundary conditions for dressing matrices of lattice systems |10j . 

From the expression ([3]), using ( }2~T|) and ( 120]) . we can write the explicit BTs as: 

_ AiA 2 (w(Xl - A|) + gfc(Aia fc - X 2 (3 k )) 
w 2 (Ai/3 fc - \ 2 a k ) 
w{\\ - Xl)a k /3 k + w 2 r k X 1 X 2 (X x a k - X 2 (3 k ) 



r k 



A?A| (Ax/? fc - X 2 a k ) 



where a k and f3 k are given by the expressions ( 12"U]) and f[2"T]) . with the boundary con- 
ditions fixed by the quadratic equations (|23|) . Before to discuss how to determine the 
sign in the solution of the quadratic equations (1231) . let us make a remark. 

A natural question can arise if one consider other type of boundary values for the 
variables a k and j3 k : one expects that the transformations will me no more energy 
preserving; we will return on this point in section [5] showing indeed how in this case it 
is possible to obtain the one-soliton solution from the "vacuum" solution (q k = r k = 
for every k). 

Now we are interested in the Hamiltonian flows generated by the transformations 
(pH|) . as outlined in section [2j We must check: 1) if the transformations are connected 
to the identity for some values of the parameters; 2) if the spectrality property holds 
[T2] . Let us see the first point. We take opposite signs in the solutions of the quadratic 
equations ( 12Q]|2~T|) for a± and fli (or a_oo and /3_oo in the case of an infinite lattice), 
that is: 



«i = w 



L 22 — Ln + A 



2L 



12 



01 = W 



-22 



L\i - A 



x=x 1 2 ^12 



A=A 2 



(25) 



with A 2 (A) = Tr(L(A)) 2 - 4det(L(A)). 
Also, let us pose: 

A 1 = e ^, X 2 = e~ tl+W , w = ie 10 (26) 

In the limit /i — > 0, from ( 1241) we obtain q k = q k , r k = r k for every k. So the 
transformations are connected to the identity. 

Let us check the spectrality property. We need to find the generating function of 
the canonical transformations f )24p . This is given by the function F(r,f) satisfying the 
equation 

dF(r, f) = V Hl ~~ qkh) dr k - ln(1 " qkTk) dr k . (27) 

V T k r k 
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Indeed, directly from the Poisson structure (fT3"|) . it is possible to check that _ ln ( 1 ^ nTn) dr n 
is a canonical one-form. From Appendix |B] the function F(r,r) is given by 

F(r, f) = ^ fkih^l) + fk{fk^\) + 9k(f k , r k+1 w 2 ) + g k {r k \{\l, r k _ 1 w 2 ) + 

- #fc(r fc f fc+ iAiA2,r fc r fc+ iW 4 ) + \n(r k w 2 ) ^ln(r fc _i^ 2 ) + ^ ln(r fc «; 2 ) + ln(r fc+ iw; 2 )^ + 
-ln(A 1 ) 2 -ln(A 2 ) 2 + J™ Hl ~ z) dz, 

(28) 

where we define the functions f k (a) and g k (a, b) by the relations 

r \n( z + r k w 2 ) . r Hz - b) , 9Q , 

fk{a) = , 9k{a,b) = . (29) 

These integrals can be explicitly evaluated in term of logarithm and dilogarithm func- 
tions. 

The variable canonically conjugate to the parameter /i is, by definition, the deriva- 
tive of the generating function ( 128]) with respect to \i. Inserting (126]) in ( 128]) and taking 
the derivative, we obtain 



dF 



2y> ■ (3°) 

This expression has to be evaluated explicitly in terms of the untilded variables r k and 
q k . As shown in Appendix [Bj the result can be written as 

s dF /Tr (L(e-» +i9 )) - A(e~^ +W )\ , s 

•^(r, q,fi,8) = -qj- — 2 In ( ^ + A( U) ) ' (31 > 

where A (A) is defined in (125]) . It is now evident that indeed this expression depends on 
the dynamical variables only as combinations of conserved quantities, since the trace 
and the determinant of L(X) are time-independent. From equation f l3~T]) . comparing 
with ([5]) and ([7]), we get the following two propositions: 

Proposition 1. The BTs fpH]) are the integral curves, with parameter /i, of the non- 
autonomous Hamiltonian system of equations governed by the Hamiltonian ^al(^, Q, 6 

dr k _ _ d$AL{r,q,fi,0) 

— = 1 - q k r k ) — , 

on dq k 

dq~k n ~ ~ ,d$ AL (r,q,fJ,,6) 

^7 = - (1 -^ } W h • 

The integral curves are identified by the initial values, given by r k \^ =Q = r k and q k \^,=o = 
Qk- 

Proposition 2. The discrete trajectories obtained by iterating the map f )24p lie, 
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for every choice of initial values (r k ,qk) and of parameters (fx, 9), on the respective 
trajectories of the following Hamiltonian system: 

dr k &H 

-m ={1 - qkrk) d q - k 



where H= / $^(r, q, X,9)dX. (33) 

Taking into account ()3ip . after some manipulations the formula for the interpolating 
Hamiltonian Ji (13"3"]) can be also written as: 

H = -2 arccosh | _5^2L ) cZA. (34) 

V-ah-w \2v/det(L(e A )) / 

A last remark: it is possible to get the reduction r k = —q* k in the formulae (124"]) 
for BTs, relevant for physical applications. In the previous formula and in the rest of 
the paper we use the asterisk to denote complex conjugation. It is easy to show that 
if r k = —ql, then the elements of the Lax matrix evaluated at A = Ai = e At+ie and 
A = A2 = e~^ +ld obey the following relations: 

-^1,1(^2) = ^2,2(^1)) ^1,2(^2) = — ^2,1(^1)1 

£2,1^2) = —L\ 2 (X\), £2,2^2) = -^i,i(Ai). 

These relations entails A 2 (Ai) = (A 2 (A2))* for the function A(A) defined in ( 125]) . The 
expressions for ol\ and f3i, or a_oo and /3_oo in the case of an infinite lattice, are explicitly 
given by ( 125]) . By choosing A(A2) = — (A(Ai))* we obtain f3\ = and, from the 

recursions ( 120H2T1) . (3 k = — -r for every k. Inserting these last equations in (I2~4"j) one has 

^k 

f k = —ql provided that r k = —q k . 

Taking into account f )26l) . the reduced BTs for the variables q k read: 



(cosh(^)(|a fc 


2 + 1) + sinh(/i)(|a fc 


2 -l))q k ~ 


h2ie 2ie sinh(2 / u)a* 


cosh(/i)( \a k 


2 + 1) - sinh(/i)( \a k 


2 -l) 



where in the definition of the variable a k ( 120]) and its boundary value (125]) the reductions 
j3 k = — -V and r k = —ql must be taken into consideration. 

^k 

5 Numerical experiments and an example 

As a first example, we take three interacting bodies. The trace of L(X) (|lip is now 
explicitly given by: 

Tr(L(A)) = A + (rig 2 + r 2 g 3 + r 3 qi)X H h -77. 

A A 13 

Taking into account also the expression for the determinant of L(X), the three integrals 
of motion are given by: 

Hi = nq 2 + r 2 q 3 + r 3 q 1: 
H 2 = qir 2 + <?2?~3 + q 3 n, 
H 3 = (1 - rigi)(l - r 2 q 2 )(l - r 3 q 3 ). 
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Figure 1: A continuous and discretised trajectory for |gi 



According to formulae ( 133]) and (134]) the discrete flow denned by the BTs ( 124]) lies 
exactly on the continuous flow generated by the Hamiltonian % given by 



The three constants of motion c±, c 2 and C3, given by formula ([8]), are explicitly 
defined by the expressions 



where A (A) is defined in (1251) . Note that in the definition of C3 there is an extra factor 
H3 in front of the partial derivative term. In this way we take into account the fact that 
in the physical Hamiltonian H% appears as the argument of a logarithm (see eq. (JHj)). 
With this choice the continuous flow discretised is given by C\Hi + C2H2 + C3 ln(i/ 3 ). 

The numerical values of the constants c±, C2 and C3 are fixed by the values of 
the initial conditions and those of the parameters /i and 9. In figure [1] we plot the 
continuous trajectory of \qi\ 2 (in gray) and the iterations of the map (124"]) (black dots), 

corresponding to the initial conditions q± = q2 = qs = d\, and = —q* k , k = 

1,2,3. In this case, the same plot could be obtained by iterating the map fl35l) . The 
values of the conserved quantities are Hi = H 2 = — | and H3 = 6; further, we chose 
the values n = 0.05 and 9 = 0.1, giving d = 0.00417 - 0.041645i, c 2 = -c\ and 
c 3 = -0.02234L 

In the next example we take = 16. As initial conditions we chose q k = 0.4 + 
0.1 cos(— n + 27r^p) and r k = —ql, corresponding, in the continuum limit, to q(x,t = 




(36) 




(37) 
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0) = 0.4 + O.lcos(x) for x G (— 7T, 7r). The values of the parameters are = 0.1 and 
9 = 0.02. In figure [2] we report the plots of \q k \ 2 for k = 1..16, obtained from the map 
(1351) . A soliton structure emerges. 

In this particular example we have sixteen conserved quantities and a BTs with two 
parameters. The constants of motion c k , given by ([8]), also depend on the parameters 
of the map. If we had sixteen parameters, by picking carefully their values it would 
be possible to obtain the particular combination of the constants c k corresponding to 
the physical flow ( fl4l) . A BTs with 2m parameters can be obtained by composing m 
times the elementary BTs ( 1351) . However, the corresponding interpolating Hamiltonian 
will no longer be given by (I34|) . To generalize equation (I34I) to the general case of BTs 
with 2m parameters appears to be difficult, but is an essential step to provide an exact 
discretisation of the physical flow (fl4l) . 

Our last example is about the "classical" use of BTs. We want to obtain a non- 
trivial flow from the vacuum solution q k = r k = 0. Considering the same boundary 
conditions on the variables a k and j3 k as discussed in the previous section, that is 
a N+1 = ai and (3 N+1 = ^ or lim^oo a k = linn^.oo a k and lim^oo f3 k = limfc^.oo f3 k , 
we get only the trivial flow q k = q k = and r k = r k = 0. This is clear from the 
fact that in this case the Hamiltonian function $ ( 13T|) is independent of the dynamical 
variables, is just a constant. By a more physical point of view it is a consequence of 
the fact that the transformations, with the aforementioned boundary conditions for a k 
and f3 k , are energy conserving. To obtain a non-trivial flow from the vacuum solution 
we must relax these conditions. In the following, for simplicity, we will retrieve the 
original notation Ai and A2 for the parameters of the BTs ( 126]) . By taking q k = r k = 0, 
the solutions of the recursions f l20|) and ff2T]) are given by a k = a X^ 2k and (3 k = PqX^ . 
In this example is no more a "time" but just a parameter. Instead, the terms «o and 
/3 will contain the time dependences. To find out the expressions for ao(t) and /?o(£) 
we must look at the Poisson structure of the model. Just as an example we can take 
the well known semi-discrete version of the zero curvature condition corresponding to 
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the equations of motion ([9]) [2], [9]: 



\ — — Tk-lA J2 — 1 _ ^kQk-l J 

where plays the role of the local Lax matrix associated to L k as in equation ([1]). 
Note that from the compatibility of equation ([[]) and equation ([2]) we can write: 

D k = M k D k - D k M k (39) 

If we evaluate ( 139|) in A = Ai, remembering that D fc (Ai)|f2^, >= 0, we have 



0. 



D k pl >= -D k \tl\ >= -D k M k pl >^D k > —M k \Vt\ > 
Since fl\, the kernel of D k , is uniquely defined up to a scalar, we obtain 

\n\ >=M k \nl>+c\nl>, 

where c is some scalar. From the previous equation we explicitly find 

c = i(A - 1), a (<) = e A i , 
where a is a constant. Analogously, for /3o(t) we have 

= e 

Inserting these expressions in (124)) and taking into account the expressions (|26|) . we 
obtain for and fv 

2sinh(2^)e i ^ (n ' i) 

Qk 



w ( e b+r,(n,t) _ e a- v (n,t)} ' _ I 0( n> t) = (3 + 2n)# + 2(1 - cosh(2/x) cos(20))t, 
2wsinh(2/i)e a+b e-^ (n ' i) W1 | r/(n,i) = (2n + 1)// + 2tsinh(2//) sin(20), 

^gb+r?(n,t) ga— r/(n,t) ' 



i.e., a one soliton solution. It can be shown that this solution coincide with the simplest 
one found through the inverse scattering method (see e.g. [2], formulae 4.5-4.6). 

A Existence of the matrix L(X) in the case of infinite 
support. 

To define the monodromy matrix for |A| = 1 we can consider the recurrence F n+1 = 
L n+ iF n , where the matrices L n are defined in fllUI) . To simplify the notation we will 
set 

/ A n (X) B n (X) 
nlJ "U„(A) D n (X) 
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The elements of the matrices F n (X) then satisfy: 



^4n+i — XA n + q n+ iC n , 

Cn+1 = ~^ + r n+lA n , 



Dn+1 — —j^ + r n+lB n , 

B n +i = XB n + q n+ iD n 



(40) 



The previous recurrences are formally the same as those defining the Jost solutions 
in the inverse scattering method (see e.g. jl], [6]): we can adopt the well known 
arguments on the existence of the Jost solution. In this appendix we will mainly follow 

m 

By denoting with G 1 the Green function for the recurrences solved by the sequences 
A n and B n , and with G 2 the Green function for the recurrences solved by the sequences 
C n and D n , we have 



G n +ik — ^G nk + 5 n ,k, 
where 5 n>k is the usual Kronecker delta function. Passing to the Fourier space 



(41) 



Ztxi 



\w\=l 



™ n 1 G i w , k dw, 5 Hjk = — 

Ztxi 



w 



n—k—l 



dw. 



\w\=l 



and inserting the above formulae in (141 p we find 

1 



G,, 



w,k 



w k (w — A) 



r 2 

^w,k 



w k (w — A x ) 



Since |A| = 1, we take a path avoiding the poles but enclosing them into the unit circle, 
obtaining: 



G„ 



n,fc 



yi-k-1 <n 

k > n 



G 2 



ri.k 



X k- n+1 

k > n 



(42) 



so that we can write: 



n-l 



A n — foX n + A n k 1 q k+ \C k , 



k=—oo 
n-l 



(43) 



fc=— oo 



where / and g are two arbitrary constants. The solution to the equations S3] is 
formally given by two Neumann series: 



An — 2^ f™ y 

p=0 

oo 

p=0 



with < 



n-l 



In — / , X Qk+ig k , 9 n - S'OA • 



k=— oo 
n-l 



(44) 



p+l _ \ ^ \k— n+1 rp fO _ f \n 



k=— oo 
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Using the result 

m / \ rn+l 



= — OO \j'=— oo / \fc=— oo / 



k=— oo V?=— oo 

valid for sequences {&&} with finite i 1 norm (see jl], Lemma A. 2), it is easy to recursively 
obtain the bounds: 

(n \P / n \P / n \ P+ 1 / n 

E n it E w > i^i^^ifo+w E 1*1 £ E w 

I* S W i ( J>|)'i ( £n)'. ^ <_ (±jr k f^± M 

From the previous expressions it follows that the Neumann series (PHI) converge for 
|A| = 1 also in the limit n — > oo. In exactly the same manner we can obtain the 
solution for the recursions on the right side of (140]) . 



B The generating function formula 

From the formula ( 1271) it follows that the generating function F(r, r) of the canonical 
transformations ( 1241) . solves the system 



dF , . 
rfc^— = -ln(l - q k r k ), 

8F < 45 > 
r k ^z~ = ln(l - q k r k ). 
dr k 



To express the right hand side of the previous equations as functions of r and f, we 
proceed in the following way. First, we solve the relation fl3]) for the variables q k , q k , 
j3 k and (3 k +i- All these four variables will be functions of r k , r k , a k , a k +i- We find the 
relations 

w{\ 1 wr k + a k - \\a k ) \\\ 2 a k (w 2 r k + \lr k ) 

qk = : , Pk 



\ia k a k+ i AiA|(w; 2 r fc + \\f k ) + w(\ 2 2 - \\)a k '' 

_ _ \ 1 (\ 1 \ 2 ,r k + w\ 2 ] a k+1 - wa k ) X 2 a k+1 (w 2 r k + \\f k ) 

Qk — 5 , Pk+l 



w 2 a k a k+l Xi(w 2 r k + A|f fc ) + w(\ 2 2 - \\)a k+1 ' 

(46) 

Substituting the index k with the index k + 1 in the expression for /3 k and comparing 
with the equivalent expression for /3 k +i, we obtain the following expression for a k+ \: 

a ^ ^ = Ai (w 4 r k r k+1 - \\\ 2 2 f k f k+1 ) 
k+1 w (Xl(Xlf k+1 - f k ) + w 2 (\\r k+l - r k )) ' 

The corresponding expressions for a k can be found by substituting the index k with 
k — 1 in ( l4"Tj) . Inserting the formulae for a k and a k+ \ in the expressions ( 146|) . we get q k 
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and q k explicitly in terms of the variables r and r. Collecting all together, the system 
(|45l) reads: 

dF_ = _ (w 2 r k + \\f k ){w 2 r k + Xlf k )(w 2 r k+1 - f k ){w 2 h-\ - AgA|r fc ) 
dr k (Xj\lf k fk-i - w*r k r k -. 1 )(\l\lr k r k+1 - w A r k r k+1 ) 

~ 9F_ = Q 2 r fc + A 2 f fc )(w 2 r fc + A%f fc )(w 2 r fc - f k _ 1 )(w 2 f k - X 2 X 2 f k+1 ) 
df k (\\\lr k r k -x - w 4 r fc r fc _i)(AfA|f fc f fc+ i - w 4 r fc r fc+ i) 

The formula (128]) for the generating function F can be checked with a direct calculation. 
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